Mon. Not. R. Astron. Soc. 000,[T]l9](2010) Printed 19 January 201 1 (MN KTbX style file v2.2) 



Star Cluster collisions - a formation scenario for the Extended 
Globular Cluster Scl-dEl GCl 



o 

(N 

1—5 



a 

Oh; 

o 



> 

oo 

m 
m 



P. Assmann^ *, M.I. Wilkinson^ t, M. Fellhauer^ ij:, R. Smith^ 

' Departamento de Astronomia, Universidad de Concepcion, Casila 160-C, Concepcion, Chile 
Department of Physics &■ Astronomy, University of Leicester, University Road, Leicester LEI 7RH, UK 



19 January 2011 



ABSTRACT 

Recent observations of the dwarf elliptical galaxy Scl-dEl (Sc22) in the Sculptor group of 
galaxies revealed an extended globular cluster (Scl-dEl GCl), which exhibits an extremely 
large core radius of about 21.2 pc. The authors of the discovery paper speculated on whether 
this object could reside in its own dark matter halo and/or if it might have formed through the 
merging of two or more star clusters. In this paper, we present N-body simulations to explore 
thoroughly this particular formation scenario. We follow the merger of two star clusters within 
dark matter haloes of a range of masses (as well as in the absence of a dark matter halo). In 
order to obtain a remnant which resembles the observed extended star cluster, we find that the 
star formation efficiency has to be quite high (around 33 per cent) and the dark matter halo, if 
present at all, has to be of very low mass, i.e. raising the mass to light ratio of the object within 
the body of the stellar distribution by at most a factor of a few. We also find that expansion 
of a single star cluster following mass loss provides another viable formation path. Finally, 
we show that future measurements of the velocity dispersion of this system may be able to 
distinguish between the various scenarios we have explored. 
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1 INTRODUCTION 

The general understanding of low luminosity stellar systems is that 
they come in two distinct flavours. On the compact side we have 
the globular clusters whose properties are consistent with a purely 
steUar composition with no need to invoke dark matt er (DM) to ex- 
plain their dynamical masses (e.g. lLaneet"ai]|2010l) . On the more 
extended side, we have the dwarf spheroidal galaxies (dSph) which 
have similar luminosities to globular clusters but are much more ex- 
tended and show enhanced velocity dispersions. These entities are 
believed to be highly DM dominated, typically havi ng estimated 
mass-to-light (M/L ) ratios of several lOs to 100s (e.g. lMatejl 19981 ; 
IWalker et al1l2009l) . with values as high as 1000 being claimed for 
some of the faintest systems (e.g. lStrigari et al.ll2008l : lKievna et al] 
These high M/L values make the dSphs the mostDM dom- 
inated objects known in the local universe (see e.g. lMated[T998l 
for a convenient review). Between th ese two re gi ons of parameter 
space there seems to be a clear gap ( iOilmore et alj2007i) in the dis- 
tribution of effective radii. This gap is visible over a wide range of 
absolute magnitudes, i.e. for several orders of magnitude in stellar 
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With the disco very of the new ultra-fa int galaxies around the 
Milky Way (see e.g. lSelokurov et al]|2010l . and references therein, 
for the latest discoveries), there seems to be a closing of the gap at 
very low stellar masses, which could be due to the tidal disruption 
of those small systems as they orbit the Milky Way. In this region, 
the distinction between star clusters and dwarf galaxies becomes 
blurred - in the absence of detailed chemical information about the 
stellar populations and star formation history, we can not be sure 
whether a particular object is a dissolving star cluster or a disrupt- 
ing dwarf galaxy. 

But this faint, low-mass end is not the only region of the 
diagram where the two distinct populations of stellar systems 
come closer to each other, reducing the gap between their half- 
light radius distributions. At the high-luminosity end, we have the 
population of ultra-compact dwarf (UCD) galaxies, which have 
masses up to 10* Mq and effective radii of 20-40 pc. Their for- 
mation mechanism is still under debate. Do we see some freak- 
ish, high-mass outl iers of the globular cluster mass-function (e.g. 
lMieskeetal]|2002l) ? Are they the product of many star clusters, 
which have form ed in a dense star forming region and then subse- 
quently merged jpellhauer & Kroup3 '2002')? Or are they instead 
the stripped nuclei of dwarf galaxies (Bekkietal. 2001)? They 
clearly fall on the star cluster side of the above-mentioned gap 
but are more extended than typical clusters, suggesting a relation 
to dwarf galaxies. Further, some of them show evidence of M/L- 
ratios which can not be explained by simple stellar population mod- 
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els ( lHaseganetZll2005l) . The jury is still out on whether or not 
UCDs have a dynamically significant DM component. 

Even in the 'normal' range of globular cluster luminosities, 
we now know of objects which are more extended than the regu- 
lar globular cluster family. Regular globular clusters have effective 
radii of order a few pc. However, there is also a population of clus- 
ters with effective radii in the range 10-20 pc. We see those ex- 
tended globular clust ers (EGC) as 'faint fuz zy' star clusters in the 
discs of SO galaxies ( Larsen & BrodidllOOOl) and as Huxor objects 
in the Halo of M3 1 jHuxor et alj2005t) . 

Extended globular clusters were first found by iHuxor et al.l 

( l2005h in M31. The authors searched for old objects in the halo of 
Andromeda and found objects which are similar to globular clus- 
ters in colour and luminosity, but with half-light radii of more than 
30 pc. Their relative dista nces to the c e ntre o f M31 are between 
15 and 35 kpc. Similarly. iHwang et al.l ( l2005h found the first ex- 
tended globular clusters in the halo of the Local Group dwarf ir- 
regular galaxy NGC 6822. The authors reported that the projected 
distances of these EGCs from the centre of NGC 6822 are close to 
13 kpc, and estimated their half-light radii to be larger than 20 pc. 
This kind of object was previously unknown and could imply the 
existence of an intermediate population of objects between usual 
GCs and dSph galaxies. Further discoverie s have been added to 
the original observations l lHuxor et al.1l2008l) . but, it is still unclear 
whether these are really a distinct population, or rather the large 
effective radius tail of a continuous distribution of cluster sizes. 

Meanwhile, EGCs have also been found in many other galax- 
ies. Extended globular star clusters (EGC) are clusters of stars for 
which the half-light radii (rh) exceed 8-10 pc. Their occurrence 
seems to be al eatory: they represent less than 3% of the clusters 
in NGC 5128 |Gomez & Woodlevll2007h . 9% of the Milky Way 
clusters I Harris 19961), an d are more than 20% of the clusters in 
M51 dChandar et alj|2004h . Even though EGCs are usually found 
at large galactocentric radii, they can also be observed in dwarf 
galaxies (Da Costa et al. 2009). This additional dichotomy within 
the star cluster population is also not well understood. As in the 
case of UCDs, a possible scenario for their formation is again the 
merging of t wo or more star cl usters, which have formed in a con- 
fined region teriins et alj|2009h . 

In this work we focus on the first EGC f ound in a dwarf ellipti- 
cal galaxy. In the Sculptor group of galaxies, lOa Costa et al. 
found an EGC (GCl) around the dwarf elliptical galaxy Scl-dEl 
(Sc22). It has a half-light radius of 21.8 pc and an absolute mag- 
nitude of Mv = -6.7. Another reported characteristic of this EGC 
is that its stellar population seems to be indistinguishable fro m the 
popul ation of the parent dwarf. By considering a king profile 
Il962h for the V-band surface brightness profile of this EGC, the 
authors were able to measure a core radius of 21.2 pc, a concen- 
tration index c= 0.65 and a central surface brightness of 23.1 V 
mag arcsec"^. Observationally, the half-light radius is defined as 
the projected radius at which the surface brightness has dropped to 
half its central value. In the theoretical literature, the term core ra- 
dius is used to refer to the natural scalelength of the model under 
consideration, for example, the King model. It is important to note 
that these definitions do not, in general coincide and in real systems 
their temporal evolution may be qualitatively different due to pro- 
cesses such as mass se gregation of stellar reimiants with hig h mass- 
to-light ratios (see e.g. lHurle\i[2007l ; IWilkinson et alj200a) . In this 
paper, we estimate the core radii of o ur models by fitting Ki ng pro- 
files to them in the same manner as in lDa Costa et al. Thus, 
our quoted core radii can be compared directly to that of Scl-dEl 
GC 1 reported in that paper. 



By means of numerical simulations we want to investigate 
whether (1) this EGC could have formed in a DM halo of its own 
(a scenario in which all globul ar clusters may re s ide iii their own 
DM halo was proposed by ,Mashchenko & Sillsl l l2005h ): (2) this 
EGC could be the product of the merging of two (or more) star 
clusters. We combine these two scenarios and simulate the merger 
of two star cluster s orbiting in a DM hal o. A similar scenario has 
been proposed by lAssmann et al. l l2O10h as a possible formation 
mechanism for the luminous components of dwarf galaxies such as 
dSphs. 

In Section|2]we describe the setup of our simulations in detail. 
In Section[3]we present and discuss our main results, namely that no 
additional DM is needed to explain the properties of GC 1 . Finally, 
in Section|4]we summarise our conclusions. 



2 SETUP 

Initially, we place two star clusters at a distance of 0.5 kpc from 
each other, at x = ±0.25 kpc. While this choice is arbitrary, it en- 
sures that the two clusters are well-separated at the beginning of our 
simulations. Each star cluster is represented by a Plummer sphere 
with a Plummer radius (effective radius) of 1 1 pc and is initially 
truncated beyond a cut-off radius of 44 pc, because less than two 
per cent of the cluster mass resides beyond this radius. The clusters 
begin the simulations in equilibrium, but we then mimic the effect 
of early gas expulsion by artificially reducing the mass of the indi- 
vidual star particles during the first crossing time of the star cluster 
(10 Myr). This mass-loss represents not only the loss of gas which 
has not formed stars, but also mass-loss from the rapid evolution of 
massive stars, although w e note that we do not use a s ophisticated 
mass-loss algorithm as in iDabringhausen et alj | |2010| ). and in the 
subsequent evolution we do not take mass-loss due to stellar evo- 
lution into account. The possible impact of mass-loss from stars on 
the late-time evolution of a merger remnant at the centre of a dark 
matter potential well is itself an interesting issue to be explored 
elsewhere. 

The initial mass of the star clusters is varied according to the 
adopted star formation efficiency (SEE) in order to obtain final 
masses (after gas expulsion) of 2.5 x lO'* MO in all simulations. 
We consider three values for the SEE, namely 10, 33 and 100 per 
cent which lead to initial masses of our star clusters of 2.5 x 10^, 
7.6 X 10"* and 2.5 x 10"* Mq, respectively. Each cluster is represented 
by 10^ particles. 

Our model clusters are placed inside spherical DM haloes of 
different masses and density profiles. As we are primarily interested 
in the early evolution of the merger remnant, we treat the halo as 
an isolated DM halo, ie.we do not simulate the dark matter halo of 
the dwarf elliptic al galaxy Scl-dEl (Sc22). We use e ither a cored 
Plummer profile l |Plummej|l91ll ; lAarseth et al.lll974) with a scale 
length of 500 pc or a cusped NEW p rofile jNavarro. Frenk & White! 
ll997l : lDehnen & McLaughliiJl2005h with a characteristic radius of 
500 pc. The mass of the halo enclosed within 500 pc (Mdm) is 
varied over many orders of magnitude namely: none, 5.6 x lO"*, 
5.6 X 10^ 5.6 X 10^ 5.6 x 10' and 5.6 x 10** Mq. The halo is mod- 
eled using 1,000,000 particles. As a cut-off radius for our halo we 
take 2.5 kpc in the Plummer cases and the virial radius for the NEW 
profiles. Due to the restrictions of fixed scale-length and fixed en- 
closed mass, the concentrations of our NEW models vary from 2. 1 
up to 89.5. In Table [T] we give an overview of the main parameters 
of the dark matter haloes we choose. 

In the case of the very low mass halo (5.6 x 10"* Mq) we also 
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Table 1. Initial parameters of the dark matter haloes used. The first column 
denotes the halo profile, which is either Plummer (P) or NFW (N). The 
second column gives the mass enclosed within 500 pc. In the third column 
we give the scale length of the halo, either the Plummer radius or the scale 
radius of the NFW profile. The fourth column gives the cut-off radius of the 
halo which is chosen to be either 5 Rpi or the virial radius of the NFW halo. 
Finally, the last column denotes the circular speed of the halo at the initial 
position of the star clusters. 
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5.52 


N 
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Table 2. Initial relative velocities of the star clusters in km s 
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V,, 


V- 
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0.00 


0.00 


0.00 


2 


0.00 


±0.10 


0.00 


3 


0.00 


±0.25 


0.00 


4 


±0.25 


0.00 


0.00 



consider different initial relative velocities for the star clusters. The 
velocities are shown in Table. |2] In the first case the clusters are 
initially at rest. In the second and third case they have transverse 
velocities, while in the last case, they have velocities directed to- 
wards each other. The choice of velocities is arbitrary, but we show 
later that our results, to first order, do not depend strongly on this 
choice of the initial velocities. 

We also perform a different set of simulations where we use 
an initial model for the two star clusters which has a scale-radius 
(Plummer-radius) of only 4 pc. This is a rnore s tandard value for 
young star clusters (see e.g. lWhitmore et al .1 19991 and many more). 

Finally, we show a suite of simulations in which we place a 
single star cluster in the centre of the halo and follow its evolution, 
without a merger, due to the expansion produced by a low SFE. 

To p erform our simulation s we use the particle-mesh code Su- 
PERBOX ( lFellhauere"tai]|200(]|) . This code has many advantages. 
One of them is that we can use an arbitrarily high number of par- 
ticles to model the cluster and the dark matter. Two-body effects 
are suppressed in a code of this type and only the smoothed poten- 
tial is taken into account. Assuming that the simulation results in a 
merged object which resides in the centre of the DM halo, thereby 
having a large relaxation-time (of order a Hubble time or longer), 
this is an appropriate approximation since we would expect that bi- 
nary formation and two-body encounter would have only a minor 
effect on our results. 

Another advantage is that Superbox allows us to set different 
levels of resolution for the different levels of high-resolution grids 
used in the code. Therefore the code provides the correct resolu- 
tion in the places where it is most needed. For our models, we use 
resolutions in the highest resolution grid of 16 pc for the DM halo 
and 0.4 pc for the star clusters. This means that the internal forces 
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Figure 1. Plots showing the properties of the remnant resulting from a clus- 
ter merger within a low mass cusped (NFW) DM halo of mass Mdm = 
5.6 X K)-' [Mq] and SFE for the star clusters of 33% (see Table|4). The sur- 
face density profile of the EGC is shown in the top panel. The (red onhne) 
crosses are the data points, and the solid curve is the King profile fit. It has a 
core radius Tc = 21.5 pc and a tidal radius ri = 196 pc. In the bottom panel 
the velocity dispersions are shown. The (blue online) circles are the line- 
of-sight velocity dispersion and the (red online) squares are the 3D- velocity 
dispersion. 



within a star cluster, the forces between the two star clusters when 
they interact and the forces from the star clusters acting on the halo 
have very high resolution, and the overall force of the background 
potential of the halo has sufficient resolution to resolve the orbits of 
the star clusters and the tidal forces acting on them, while keeping 
computational costs low. 

In this code, two-body effects, like binary formation and two- 
body relaxation are neglected, since they do not play an important 
role if the final object resides in a DM potential. 

3 RESULTS 

3.1 Cluster merger without DM 

The first simulations we perform are of two star clusters which 
merge under their own gravity without a DM halo present. In Ta- 
ble[3]we show the properties of the merger object after 10 Gyr of 
evolution. This long time-span is chosen to ensure that the object 
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Table 3. Properties of the merger objects in the no-DM simulations. The 
first column gives the SFE. The second, third and fourth columns show 
the parameters of a King profile fit to the resulting merger object, namely 
the core radius r^, the tidal radius and the central surface brightness //q 
assuming a stellar mass to light ratio of 1 .0. This choice is arbitrary, because 
the real value is u nknown, but i n agreement with visual M/L ratios of many 
old globulars (e.g. lHarriJl99^ . 
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has h ad time enou gh to virialize and evolve slowly. We fit King pro- 
files ( lKing|ll962h to the final merger objects for this and all mod- 
els studied because these profiles are widely used in observational 
astronomy and, in particular, were used in the discovery paper of 
EGC Scl-dE GCl jPa Costa et"ai]|2009h . In reality, we would not 
expect a merger remnant to exhibit a profile with a sharp tidal trun- 
cation radius, as in a King profile. Therefore the nominal value of 
the King tidal radius should be regarded more as a third parameter 
of the fit rather than a physical property of the remnant. 

In the SFE= 100 % case, the fit gives a core radius = 1 pc. 
This is obviously a smaller value than the initial Plummer radius but 
this does not imply that the final merger object is denser than its two 
constituents. Fitting a Plummer profile to the merger object gives 
a fitted Plummer radius of about 1 1 pc. So the apparent reduction 
in only arises due to the different physical meaning of the scale 
radius in the two profiles (King and Plummer). 

If we introduce super- virial, expanding systems by increasing 
the initial virial mass of the star cluster (which then gets artificially 
reduced over one crossing-time), thereby mimicking a lower SFE 
effciency of only 33 % (or 10 % respectively) and the expansion due 
to the gas-expulsion, we do not find a merger object at all. The two 
star clusters simply disperse and build a very extended distribution 
of stars, which in the presence of an external host galaxy would 
have no chance to survive. 



Table 4. Results of the simulations of cluster mergers in low mass haloes. 
In all cases, the mass of the halo (enclosed within 500 pc) is 5.6x 10* [Mq]. 
The first three columns indicate the profile of the DM halo (P for a Plummer 
profile, N for a NEW profile), the choice of initial cluster relative velocities 
(with notation according to Table|2), and the adopted SEE. Columns four to 
six show the parameters of the fitted King profile to the final merger object 
(i.e. core radius r^, half-light radius rj, (for the fitted King profile), tidal ra- 
dius ft and cenh'al surface brightness jiio)> while the seventh column gives 
the ratio of DM to luminous matter within the central 50 pc of the object. 
This can be transformed into a regular M/L ratio by adding the masses of 
both components and adopting a stellar M/L ratio for the luminous compo- 
nent. For example, having a mass-to-mass ratio of unity and adopting a M/L 
for the luminous component of unity as well, we get a total M/L of 2. 
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3.2 Cluster merger with DIVI 

We now place both star clusters in a common, very low-mass DIM 
halo (Mdm = 5.6 x lO"* IVIo) and test for the influence of the en- 
counter geometry. The results of these simulations are given in Ta- 
ble O For the SFE= 100 % cases, the core radii of the resulting 
EGCs are independent of the initial relative velocities of the clus- 
ters, with Plummer profiles yielding EGCs with = 7 pc while 
NEW haloes lead to more compact EGCs with = 3.5 pc. We find 
two simulations, using NEW haloes, in which the SCs do not merge 
and remain as two independent structures. These are the cases in 
which we give the star clusters transverse velocities of the order of 
the halo circular velocity. In those cases the star clusters meet ofF- 
centre but are compact enough not to merge. This only happens in 
the NEW models because here the cusp ensures that the orbits of 
both clusters never lead to close encounters between the clusters. 
Because of the low mass of the halo and the low densities in the 
outer regions of the two clusters (we use steep Plummer profiles 
as inital models), dynamical friction is not strong enough to bring 
both clusters close enough to merge. However, EGCs always form 
if the SFE is 33 %. Two of these simulations, which have a Plum- 
mer DIVI halo, result in EGCs with = 22 pc. The EGC formed 



inside a NEW D]VI halo which most closely resembles the proper- 
ties of Scl-dE GCl is formed when the relative velocities have a 
small transversal component (case 2 in Table O, i.e. the two star 
clusters do not merge in a head-on collision. In the top panel of 
Figure[T]we show the surface density profile of this cluster. 

It is interesting to note that in all cases with a low SFE of 10 %, 
we do not form a merger object at all but the two clusters instead 
dissolve more quickly than they are able to merge. The result is 
a system with a diffuse lig ht distribution, more s imilar to a dSph 
galaxy than a star cluster. In lAssmann et al. U2OIOI) and Assmann et 
al. (in prep.), we explore this scenario for the formation of dSphs 
in detail. 

In Figure [21 we show in the top-left panel the resulting core 
radii as function of initial cluster relative velocity (see Table [2J. It 
shows clearly that the size of the merger object depends mainly on 
the adopted SFE and not on the geometry of the first encounter, as 
long as the two clusters are able to merge. We find no significant 
general trend in our simulations. As a minimal, second order trend 
we see, that if the star clusters meet slightly off-centre (velocity 
cases 2+3) the resulting core radii are mostly slightly larger and if 
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Figure 2. Top-left: Core radius of our merger objects as function of the adopted relative velocity case (see Table |2)- The halo mass was 5.6 X lO'* Mo m 
all cases. Triangles (blue online) show Plummer profiles and diamonds (green online) show NFW profiles for the halo with a SFE of 100 % for the two star 
clusters. Stars (red onhne) and points (magenta online) show Plummer and NFW shapes for the halo, respectively, if the SFE is 33 %. Top-right: The core 
radii of our merger objects as function of the mass of the halo. Triangles (blue onhne) show Plummer haloes and star clusters with 100 % SFE. Diamond (green 
onhne) the same SFE but with NFW haloes. Stars (red online) show Plummer haloes using 33 % for the SFE and points (magenta online) the con'esponding 
NFW cases. The results for the simulations without DM are for compaiison shown at a halo mass of 1000 because of the logaiithmic scale. The only resulting 
merger object at a SFE of 10 inside a NFW halo is shown as the square (black onhne). Bottom-left: Core radius versus halo mass for our models with compact 
star clusters. Triangles (blue) and Diamonds (green) show star clusters with 100 % SFE in haloes with Plummer and NFW profiles, respectively. Stars (red 
online) show Plummer type haloes and a star cluster with a SFE of 33 %. Circles (magenta) show the same for haloes with NFW profiles. Bottom-riglit: Core 
radius versus halo mass for our single star cluster models. Stars (red online) show Plummer type haloes and a star cluster with a SFE of 33 %. Points (magenta) 
show the same for haloes with NFW profiles. Square (black) and Triangles down (orange) show star clusters with 10 % SFE in haloes with Plummer and NFW 
profile, respectively. Black asterisk show the resulted star cluster from initial compact star cluster evolved without dark matter halo. See text for a detailed 
discussion. 



we increase the speed of the head-on encounter (case 4) we get in 
most of the cases a smaller core radius. 

As a next step, we increase the mass of the halo by factors of 
10. The results of our simulations are given in Table |5] Again if 
the SFE is 100 % we only obtain merger objects with core radii of 
only = 3-7 pc independent of the mass of the halo. As soon as 
the SFE is 33 %, it is possible to obtain an EGC which resembles 
Scl-dEl GCl as long as the halo does not become too massive, 
i.e. the DM content within the final merger object is well below its 
luminous mass. However, it is interesting to note that in the two 
cases with = 5.6 x 10^ Mq the EGCs have a concentration 



(c = log 10(r,/rc)) that is double the estimated concentrati on of 
0.65 for EGC Scl-dE GCl reported by l lDa Costa et alj2009l) . 

If we consider a mass for the DM halo of = 5.6 x 10* Mq 
and a SFE of 33%, the simulated EGCs have core radii which are 
too large (r^ = 26 - 39 pc). In other words when the DM content 
within the merger object is of the same order as that of the luminous 
matter (mass ratio of 1 .7) or significantly higher (mass ratio of 24.4) 
we get merger objects which are too extended. However, in those 
cases it may be possible to find a matching model by tuning the SFE 
to higher values. In almost all simulations with a SFE of 10 %, the 
SCs dissolve before they merge. 
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Table 5. Results of our simulations for dilferent halo masses (enclosed 
within a radius of 500 pc). The first column gives the shape of the halo 
(- for none, P for Plummer and N for NFW), the second column the mass 
of the halo and the third the SFE. The fourth, fifth and sixth columns give 
the parameters of a King-profile fitted to the data, namely the core radius, 
tidal radius and central surface brightness. The last column gives the mass 
ratio between the DM and the luminous matter within the central 50 pc. The 
results shown are mean values from up to four different realisations. 



profile 
P/N 



Ml, 
[Mo 



SFE 
[%] 



[pc] 



[pc] 



[pc] 



[mag/ 
arcsec"! 



M.,iar 





0.0 


100 


7.0 


18.5 


50.1 


21.3 


0.0 




0.0 


33 






dissolved 






p 


5.6x10-* 


100 


7.0 


18.4 


49.7 


21.5 


0.025 


N 


5.6x10'' 


100 


3.5 


11.1 


31.0 


29.4 


0.01 


P 


5.6x10* 


33 


18.3 


49.6 


135.0 


24.1 


0.01 


N 


5.6x10'* 


33 


23.1 


54.0 


143.0 


24.6 


0.06 


P 


5.6x10* 


10 






dissolved 






N 


5.6x10* 


10 






dissolved 






P 


5.6x10'^ 


100 


3.0 


11.0 


31.5 


19.6 


0.1 


N 


5.6x10"' 


100 


3.5 


12.3 


35.2 


19.6 


0.13 


P 


5.6x10^ 


33 


21.3 


85.6 


249.0 


24.4 


0.31 


N 


5.6x10' 


33 


19.3 


103.7 


313.0 


24.7 


0.06 


P 


5.6x10"' 


10 






dissolved 






N 


5.6x10"^ 


10 






dissolved 






P 


5.6x10* 


100 


4.1 


9.5 


25.2 


21.1 


1.06 


N 


5.6x10* 


100 


4.6 


29.6 


91.0 


20.8 


3.37 


P 


5.6x10* 


33 


26.8 


45.4 


112.0 


23.6 


1.69 


N 


5.6x10* 


33 


36.8 


596 


1950 


25.8 


24.4 


P 


5.6x10* 


10 






dissolved 






N 


5.6x10* 


10 


106 


198 


500 


26.0 


40.8 



P/N 


5.6x10'' 


all simulations dissolve 


P/N 


5.6x10** 


all simulations dissolve 



Simulations using a higlier mass for the DIVI halo (5.6 x 
10'' [IVIq] and 5.6 x 10* [IVIq] enclosed within 500 pc), lead to 
completely dissolved star clusters independent of the SFE used. 
In none of these simulations did we actually see the merging of 
the two clusters. This can be understood using the theory devel- 
oped in iFellhauer etal] | |2009|) . If the strength of the bacltground 
potential becomes too high, i.e. the star clusters meet each other 
with relative velocities much higher than their internal dispersions 
it results in destructive high-speed encounters rather than in low 
velocity merging. 

In the top-right panel of Figur42]we plot the core radius of our 
merger objects as function of the halo mass. We clearly see that the 
mass of the halo has only a small influence on the final core radius 
of the resulting merger object, except for the fact that we have no 
merger object in very high mass haloes. The main parameter which 
governs the results is clearly the SFE of the two star clusters. This 
SFE has to be of order 33 % to obtain a merger object which re- 
sem bles the EGC. 

iDa Costa et all l l2009h discussed whether the EGC Scl-dE 
GCl could be a dark matter dominated object. Our most impor- 
tant result is that we only get objects which resemble this EGC 
if the mass of the DIM halo within the final stellar distribution is 
negligible. Higher halo masses lead either to objects which are too 
extended (i.e. we have to use higher SFEs to counteract external 
destruction) or to no merged objects at all. To understand this be- 



Table 6. The central line-of-sight velocity dispersions of our best-fitting 
models, i.e. whose is close to that of Scl-dE GCl. All cases shown here 
have a SFE of 33 %. For a halo mass of 5.6 X 10* M© column three gives the 
velocity case, while for halo masses of5.6xl0'Mo we give the realisation 
number 





profile 


velocity 




0"L0S,mean 


[Mo] 


P/N 


case 


pc 


km S"' 


5.6 X 10* 


N 


2 


21.5 


0.78 


5.6 X 10* 


P 


3 


22 


0.94 


5.6 X 10* 


P 


4 


22 


1.17 


5.6 X 10* 


N 


4 


21 


1.68 


5.6 X 10^ 


P 


1 


21.4 


1.63 


5.6 X 10' 


P 


2 


23 


1.66 


5.6 X 10^ 


N 


2 


19.5 


1.66 


5.6 X 10^ 


P 


3 


19.6 


1.65 


5.6 X 10^ 


N 


3 


19.5 


1.62 


5.6 X lO' 


P 


4 


21 


1.64 


5.6 X 10^ 


N 


4 


20 


1.61 



havior we need to consider all the competing processes which drive 
the evolution of the merging clusters within the DIVI halo. First, it 
is important to observe that the expansion of the star clusters is 
governe d by the SFE, such tha t lower SFEs give more extended 
objects teaumgardt et alJboOTi ). Secondly, it must be emphasized 
that the influences of the DIVI halo can be divided into two distinct 
processes, that depend on where the SCs are initially located. One 
process is when the star clusters have not yet merged in the cen- 
tre of the DIVI halo. In this case, the tidal forces generated by DIVI 
halo act as an additional destructive force on the clusters, reduc- 
ing their masses and sizes. The more massive the DIVI halo, the 
more dominant is this disruptive process. Additionally, a higher 
halo mass increases the final relative velocities with which the SCs 
meet and might, therefore, inhibit the merging. This first halo pro- 
cess tends to lead to remnants which are more extended and diffuse 
than observed EGCs. The second halo process takes place when the 
merged object forms and remains at rest in the centre of the DIVI 
halo. In this case, the potential well of the D]VI halo contributes 
to the gravitational field in which the cluster stars move. As a re- 
sult, objects which might otherwise dissolve are bound by the DIVI 
halo, and the additional gravitational field makes the remnants more 
compact than they would be without the presence of DIVI. 

Our simulations suggest that mass-to-light ratios of only a few 
would be expected if the object is formed by the merger of two star 
clusters, and thus, they are consistent with there being no dynami- 
cally significant amount of DIVI within the stellar distribution of the 
EGC. However, our results also show that we do require some DIVI 
to ensure that the clusters merge to produce an EGC-like remnant, 
so that DM is an essential part of the merger scenario per se despite 
not being a significant component of the merger remnant. 

In addition to comparing the core radius and concentration 
of our merger remnants with those of EGC Scl-dE GCl, we can 
also use our simulations to look at the velocity space of the EGC. 
Table |6] shows the mean, central line-of-sight velocity dispersion 
(o"Los,mean) for the simulations which resemble the EGC. We av- 
erage the central dispersion measured along all three coordinates 
axes (the two coordinates which span the plane of the interac- 
tion and the one perpendicular to that) because the orientation of 
the sight-line towards GCl relative to any proposed merger is, of 
course, unknown. In our models, (TLos.mcaii varies from 0.78 km.5"' 
to 1.68 km .V"'. In Figure [T] we show the dispersion profile for one 
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Table 7. Results of our simulations for star clusters which are initially more 
concentrated, i.e. have a Plummer radius of 4 pc. The columns are the same 
as in Table [5] 



profile 
P/N 
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[pc] 



Hi 

[pc] 



[pc] 



/JO 
[mag/ 
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1 fin 






did not mcr^ 


5C 














did not mGr^ 


5C 




p 


5.6 X 10" 


100 


2 


25.1 


81 


19.2 


0.03 


N 


5.6 X 10" 


100 


2.8 


12.5 


37 


19.1 


0.01 


P 


5.6 X lO" 


33 


6.8 


12.1 


30.3 


22.05 


0.06 


N 


5.6 X lO" 


33 


7.2 


24.1 


68.2 


22.0 


0.02 


P 


5.6 X 10" 


10 






dissolved 






N 


5.6 X lO" 


10 






dissolved 






P 


5.6 X 10^ 


100 


2.9 


9.9 


28 


19.3 


0.13 


N 


5.6 X 10^ 


100 


2.8 


14.0 


42 


19.4 


0.2 


P 


5.6 X 10^ 


33 


12.5 




1.2e+5 


24.2 


0.35 


N 


5.6 X 10^ 


33 


10.7 




1.2e+3 


23.9 


0.7 


P 


5.6 X 10^ 


10 






did not merg 






N 


5.6 X 10' 


10 






dissolved 






P 


5.6 X 10** 


100 


4.1 


10.4 


28 


21.88 


1.1 


N 


5.6 X lO* 


100 


3.2 


33.2 


106 


20.33 


3.34 


P 


5.6 X 10* 


33 


18.7 




3.1e+4 


24.5 


0.35 


N 


5.6 X 10'^ 


33 


57.2 




1.2e+5 


26.1 


24.1 


P 


5.6 X 10* 


10 






dissolved 






N 


5.6 X 10* 


10 






dissolved 







of our merger objects. It is notable that it is possible to increase 
the central dispersion of our merger object up to almost 1.7 Ian s"' 
without having significant amounts of DIVI within the object. 

The importance of studying the velocity space of the simulated 
clusters is that the velocity dispersions can be used to distinguish 
between our various scenarios for the formation of EGC Scl-dE 
GCl. Moreover, the velocities obtained in our simulations can be 
compared with future observations to infer the presence or absence 
of dark matter in the system. 



3.3 Compact initial clusters 

For comparison, we also perform simulations with a more standard 
model for the initial star cluster, i.e. with a Plummer radius of just 
4 pc. With this suite of simulations we want to investigate whether 
it is possible to obtain EGCs like GCl without starting with star 
clusters that are already extended. 

The results in Table [7] and the bottom-left panel of Figure (2] 
show clearly that we do not end up with a merger object that resem- 
bles the properties of GCl, if the initial star clusters are too concen- 
trated. A SFE which guarantees the survival of the two clusters al- 
ways leads to merged objects which are too concentrated to resem- 
ble the EGC in Sculptor. All our merged objects show core radii of 
less than or about 10 pc. This behavior is explained by the force per- 
formed by the DIVI halo as soon as the merger SCs are settled in the 
centre. This force keeps the merged object together and more com- 
pact. The exceptions are the simulations with a 5.6 x 10* IVIo halo. 
Here again we see a completely dissolved object but this time it is 
possible to determine a centre of density. However, even though the 
fitted 'core radii' seem to match the properties of the EGC, one sees 
from the values of the 'tidal radii' that these are not really EGCs but 
rather fluffy, low-density distributions of stars. Taking these results 



Table 8. Results of our simulations with a single star cluster expanding 
because of mass-loss due to gas-expulsion. The columns are the same as in 
Table |5] For the non-DM simulations we give two simulations which have 
a SFE of 33 %. The first simulation started with a more compact Plummer 
sphere and has a Plummer radius of 4 pc, while the second one shows the 
standard 1 1 pc model. 
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[mag/ 
arcsec^] 






0.0 
0.0 
0.0 


33 
33 
10 


11.3 
59.9 


21.1 
62.7 


53.4 
139 

dissolved 


24 
27.6 




P 
N 
P 
N 


5.6 X 10" 
5.6 X 10" 
5.6 X 10" 
5.6 X 10" 


33 
33 
10 
10 


32.9 
23.9 


34.5 
42.3 


76.5 
105.6 
dissolved 
dissolved 


23.6 
24.2 


0.23 

0.03 


P 
N 
P 
N 


5.6 X 10' 
5.6 X 10^ 
5.6 X 10' 
5.6 X 10' 


33 
33 
10 
10 


22.9 
23 
47.1 
71.7 


39.2 
34.3 
84.0 
67.7 


97.2 
82.3 
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147 


23.7 
23.3 
26.3 
25.7 


0.21 
0.35 
0.55 
0.5 


P 
N 
P 
N 


5.6 X 10* 
5.6 X 10* 
5.6 X 10* 
5.6 X 10* 


33 
33 
10 
10 


20.1 
16.8 
47 
25.8 


37.9 
29.1 
51.4 
86.0 


96 
72.2 
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243 


23.1 
22.5 
24.1 
24.1 


0.8 
4.1 
1.34 
4.2 


P 
N 
P 
N 


5.6 X 10^ 
5.6 X 10^ 
5.6 X 10' 
5.6 X 10' 


33 
33 
10 
10 


17.7 
9.2 
33.9 
32.8 


29.6 
36.3 


73 
3.5e+4 

81 
4.7e+5 


22.5 
23.8 
23.01 
25.34 


5.89 
74.8 
7.78 
91.2 



into account we can rule out that the EGC has formed out of two 
concentrated star clusters. 



3.4 Abandoning the merger scenario 

To complete our survey of possible formation scenarios we now 
abandon the hypothesis that Scl-dE I GCl formed via a merger of 
multiple clusters and explore whether it is possible to expand a sin- 
gle young star cluster sufficiently due to gas-expulsion to resemble 
the observed EGC. 

Our simulations (results are shown in Table [8j show clearly 
that it is possible to get an extended star cluster just by expansion 
due to gas-expulsion alone. In our simulations we do not include 
any tidal field of the parent galaxy and therefore our star cluster 
can expand freely without fearing destruction due to external tidal 
forces. In the case without DM, we use only our standard SFE of 
33%, which has been shown in many studies (e.g. lBaumgardt et al.l 
|2007|) to be the lowest limit for a cluster to survive gas-expulsion. 
Therefore it is no wonder that the cluster expands significantly. 
However, by using higher SFE and/or initially more compact clus- 
ters it is likely that we could find a suitable model to reproduce the 
observed data. We performed one simulation which started from a 
compact star cluster with a Plummer radius of only 4 pc and was 
transformed into an EGC with a scale length of 11.3 pc, while a 
cluster with an initial Plummer radius of 1 1 pc was transformed 
into an EGC with a scalelength of ~ 60 pc. This clearly shows that 
we would be able to get a model which reproduces the observed 
data by using an intermediate initial scale-length. Applying Oc- 
cam's Razor, this might be the simplest and most plausible way to 
explain the origin of the EGC in Sculptor. However, the intention 
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Table 9. The central line-of-sight velocity dispersion of our best-fitting, 
single-cluster models, i.e. whose is close to the one of Scl-dE GCl. All 
cases shown here have a SFE of 33 %. In the first two lines we give the re- 
sults obtained in our non-DM simulations. The values of in these models 
bracket that of the observed EGC - as discussed in the text, we expect that 
a good model of the observed data could be obtained for an intermediate 
value of the initial cluster Plummer radius. 
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0.18 


5.6 X 10'' 
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1.06 


5.6 X 10' 


N 


23 


1.73 
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P 


22.9 


1.66 


5.6 X lO* 


P 


20.1 


2.27 



of this paper is mainly to test and investigate the alternative, and 
more complicated, formation scenarios which have been proposed 
in the literature. 

If we place the single cluster in the centre of its own DM halo, 
we first see that the more massive the halo gets the more likely it 
is for the star cluster to survive even for SFEs as low as 10 %. The 
higher the halo mass the more gravitational potential there is to halt 
the expansion and therefore the more compact our resulting object 
is. This trend is visible in the bottom-right panel of Figure[2] With 
our initial cluster model and values for the SFE we get the best 
matches for intermediate mass haloes of 5.6 x lO"* to 5.6 x 10* Mq 
for both kinds of halo profile. Again we see no significant trend 
with halo profile except that we expect the NFW haloes to have 
more mass in the region of interest than the equivalent Plummer 
model. However, from the general behavior of our results we expect 
to be able to find a suitable solution for any halo mass. If we had 
extended our parameter space to include even more massive haloes 
then we would likely have obtained a match even for SFE of 10 %. 
At lower halo masses we could use either initially more compact 
clusters or assume higher SFEs. 

Finally in Table|9]we show the LOS velocity dispersion of our 
best-fitting, single-cluster models. As expected, the central velocity 
dispersion rises with DM halo mass, from below one km s"' for the 
non-DM case up to more than two kms ' for the DM dominated 
one, and therefore will be the crucial parameter to decide whether 
or not this very extended object has its own DM halo, once spec- 
troscopic data are available. 



4 CONCLUSIONS 

In this paper, we explored several of the proposed formation sce- 
narios for the recently-discovered, extended globular cluster GCl 
associated with the dwarf elliptical galaxy Scl-dEl in the Sculp- 
tor group of galaxies. Even though extended globular clusters have 
been found in the haloes of many galaxies, this is the first such 
object found in orbit around a dwarf galaxy and exhibiting a very 
large core radius. 

There are two obvious scenarios for the formation of an ex- 
tended stellar system starting from one, or several, individual con- 
stituents. The first, and probably more natural, scenario is via the 
expansion induced by either gas expulsion or stellar mass loss dur- 
ing the early evolution of a normal star cluster. The other viable 
scenario is that extended objects form via the merger of smaller 



systems. A third, more recent proposal for the formation of globu- 
lar clusters is that they might reside in their own DM halo, thereby 
helping to solve the missing satellite problem. We have extended 
those speculations to the object of our study. In this paper, we have 
combined aspects of these three scenarios in order to investigate 
thoroughly which initial conditions can lead to an object which re- 
sembles the observations. 

Taking all our merger simulations into account, we see the 
following trends: 

(i) As expected, for an assumed SFE of 100%, mergers of two 
clusters without the presence of DM produce a remnant which is of 
similar density to the original clusters, while lower SFEs lead to the 
dissolution of both initial clusters without merging. Finely tuned 
initial conditions would be required to produce a merger remnant 
via this route which matches Scl-dEl GCl in terms of size. 

(ii) Mergers of two virialised clusters within a DM halo do not, 
in general, lead to a remnant which is more extended than the orig- 
inal clusters. This result is independent of the assumed halo mass. 
We therefore need to invoke additional expansion either due to a 
SFE which is lower than 100 % or due to mass loss from early 
stellar evolution to produce extended remnants. 

(iii) For high halo masses, we find that low SFEs can lead to 
the dispersal of the star clusters before they have time to merge. In 
those cases we see an extended luminous component which more 
closely resembles a dSph galaxy than a globular cluster. This co- 
in cides with the format ion theory for dSph galaxies we discuss 
in lAssmann et al. and Assmann et al. (in prep.). 

(iv) In all the simulations which led to an extended remnant 
whose appearance matched that of Scl-dEl GCl, the mass of DM 
enclosed within the volume probed by the stars of the EGC was not 
sufficient to raise the mass to light ratio significantly above the stel- 
lar M/L. Therefore, even though our simulations can not rule out 
that this extended globular cluster resides in its own DM halo we 
do not expect it to be a highly DM dominated object at the present 
day. 

(v) We have shown that our models can be tested by future ob- 
servations. In our remnants, the higher the DM content, the higher 
is the central velocity dispersion of our merger object ranging from 
0.7 to 1.7 kms"'. The next generation telescopes, for example the 
Extremely Large Telescope, may be able to measure such differ- 
ences either from the integrated spectra of the clusters or from ob- 
servations of individual cluster stars, which can then be compared 
with the predictions of the models we have simulated. 

(vi) If we merge clusters which are initially as compact as typ- 
ical globular clusters, we find that it is very difficult to form an 
object resembling the observed EGC. We need to lower the SFE 
below the general survival limit while still requiring that the clus- 
ters merge in the centre of the halo before they dissolve completely. 
Such chance m ergers are possible and hav e been found in other 
studies as well jFellhauer & Kroupall2005l) . but are unlikely to be 
a dominant formation avenue. 

Abandoning the merger scenario, we have shown that we can 
also form an object like GC 1 via the expansion due to gas-expulsion 
during the formation of the cluster. In these models, we start with 
an object forming in the very centre of a DM halo and show that gas 
expulsion produces remnant clusters whose extent depends on the 
mass of the halo - the higher the halo mass the more compact is the 
remnant. With high halo masses it is possible to get surviving and 
matching objects by using SFEs which would lead to the complete 
dissolution of the cluster in the absence of DM. Starting out with a 
more compact object leads to a final remnant which is more com- 
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pact. Thus, for any given halo mass, it is possible to match the data 
either by changing the SFE or by changing the initial size of the 
object. Again, as in the merger scenario, the observable difference 
between the models which resemble Scl-dEl GCl is their veloc- 
ity dispersion. The higher halo mass, and therefore the DM content 
within the EGC, the higher its velocity dispersion will be. 

The velocity dispersions of globular clusters containing 
Intermediate -Mass Black Holes (IMB H) at their centres have been 
discussed bv lBaumgardt et al] | |2005l) using numerical models. The 
velocity dispersion for a star cluster of similar stellar mass to our 
EGC is raised by only 0.1-0.5 km in the inner 10 per cent of 
the half-mass radius. The same would also be true for the stronger 
effect of a binary IMBH. In both cases, the increase in velocity dis- 
persion is seen only in the core the star cluster. Other numerical 
simulations show, that ebven if t here is a high binary black hole 
fraction within a globular cluster jMackev et alj|2008f) . the central 
line of sight velocity dispersion is about 1.0-1.5 km s^^ about dou- 
ble the value, but shows the regular fall-off at the outer radii as 
well. 

To distinguish between the presence of a black hole and the 
presence of dark matter distributed within the cluster would require 
the measurement of the velocity dispersion as a function of radius 
within the cluster. If the EGCs contain dark matter, then their ve- 
locity dispersions will remain higher than expected at large radii, 
not solely at the centre. For some EGcs, such observations may be 
feasible with the next generation of extremely large telescopes. 

Finally, we note that it is possible to obtain an object like Scl- 
dEl GCl just by using expansion due to mass-loss (gas expulsion 
and/or stellar winds) within a single young cluster. Fine-tuning the 
initial concentration of the cluster and, in our case, the SFE leads to 
an object which matches the data. In the light of Occam's razor this 
might be the most simplest and straightforward solution to explain 
the properties of Scl-dEl GCl. The key test of this model would be 
the survivability of the remnant in the tidal field of the host galaxy. 
This would require further simulations of the evolution of the clus- 
ter in the external tidal field following the period of expansion and 
is beyond the scope the present paper. 

We conclude that the unusually large size of the extended star 
cluster Scl-dEl GCl does not require the presence of dynamically 
significant DM in this system. We have shown that there are several 
formation paths for this object involving mergers of star clusters 
in low-mass dark matter haloes which result in a remnant that 
contains very little DM interior to the stellar distribution. It is also 
possible to form the object through the expansion of a single star 
cluster as a result of mass loss. Thus, we have demonstrated that 
the formation scenarios proposed for this system in the literature 
are plausible. To make further progress requires the measurement 
of a precise velocity dispersion for the object, as this would enable 
us to distinguish between formation models which require DM and 
those which require only baryonic matter. To distinguish between 
scenarios which involve merging of two or more initial objects and 
models which only expand due to gas-expulsion alone, one might 
also be able to look for irregular chemical abundances in the stars 
of the cluster which could point to an origin in different clusters. 
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